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Abstract. Azimuthal asymmetries in exclusive elcctroproduction of real photons are 
measured for the first time with respect to transverse target polarisation, providing 
new constraints on Generalized Parton Distributions. From the same data set on a 
hydrogen target, new results for the beam-charge asymmetry are also extracted with 
better precision than those previously reported. By comparing model calculations 
with measured asymmetries attributed to the interference between the deeply virtual 
Compton scattering and Bcthc-Heitler processes, a model-dependent constraint is 
obtained on the total angular momenta carried by up and down quarks in the nucleon. 



1. Introduction 

The partonic structure of the nucleon has traditionally been described in terms of 
Parton Distribution Functions (PDFs) of the parton's longitudinal momentum as a 
fraction of the nucleon's momentum in a frame in which the nucleon is moving at 
almost the velocity of light. These functions appear in the theoretical description of, 
e.g., Deep- Inelastic Scattering (DIS). However, in the context of the rapid theoretical 
developments of the last decade, PDFs have been conceptually subsumed within the 
broader framework of Generalized Parton Distributions (GPDs), which also describe 
elastic form factors and amplitudes for hard-exclusive reactions leaving the target 
nucleon intact [1, 2, 3]. Most often discussed are the four twist-2 quark- chirality 
conserving quark GPDs: the polarisation-independent distributions H q and E q and 
the polarisation-dependent distributions H q and E q . The GPDs H q and H q conserve 
nucleon helicity, while E q and E q are associated with a helicity flip of the nucleon. GPDs 
depend on the kinematic variables x and £, which represent respectively the average and 
difference of the longitudinal momentum fractions of the probed parton in the initial and 
final states. The variable £ is typically nonzero in hard-exclusive reactions. GPDs also 
depend on the squared four-momentum transfer t = {p — p') 2 to the nucleon, with p (p') 
the four-momentum of the nucleon in the initial (final) state. PDFs and nucleon elastic 
form factors appear as kinematic limits (t —>■ 0) and x-moments of GPDs, respectively. 
Strong interest in the formalism of GPDs and their experimental constraint has emerged 
after moments of certain GPDs were found to relate directly to the total (including 
orbital) angular momenta carried by partons in the nucleon, via the Ji relation [2]: 



This finding offers for the first time a path towards solving the 'nucleon spin puzzle' of 
how the helicities and orbital angular momenta of quarks and gluons combine to form 
the spin of the nucleon. More recent discussions have focused on the potential of GPDs 
as multi-dimensional representations of hadrons at the partonic level, correlating the 
longitudinal momentum fraction with transverse spatial coordinates [4, 5, 6, 7, 8]. 



PACS numbers: 13.60.-r 24.85. +p 13.60.Fz 14.20.Dh 




(1) 




(a) (b) 

Figure 1. Leading-order diagrams for (a) deeply virtual Compton scattering 
(handbag diagram) and (b) Bethe-Heitler processes. 

2. Deeply virtual Compton scattering 

Generalized Parton Distributions are accessible through exclusive processes that involve 
at least two hard vertices, yet leave the target nucleon intact. Among all presently 
practical hard exclusive probes, the Deeply Virtual Compton Scattering (DVCS) 
process, i.e., the hard exclusive leptoproduction of a real photon (e.g., 7* N — > jN'), 
appears to have the most reliable interpretation in terms of GPDs. In electroproduction, 
direct access to the DVCS amplitude T DV cs is provided by the interference between the 
DVCS and Bethe-Heitler (BH) processes, in which the photon is radiated from a quark 
and from the lepton, respectively (see Fig. 1). Since these processes are intrinsically 
indistinguishable, the cross section is proportional to the squared photon-production 
amplitude written as 

T| 2 I rp |2 . Irp |2 . rp rp* , rp* rp /Q\ 

I — I t DVCS I + I -L BH I + 1 DVCS 1 BH + 1 DVCS 1 BH , {*) 

I 

where T denotes the interference term. The BH amplitude Tbh is precisely 
calculable from measured elastic form factors of the nucleon and provides the dominant 
contribution in Eq. 2 in the kinematic conditions of the present measurement. These 
amplitudes depend on Q 2 = —q 2 with q = k — k! and k (k') the four-momentum of the 
lepton in the initial (final) state, the variable xb = Q 2 / '(2Mv) with v — p ■ q/M and 
M the nucleon mass, and t. In addition, the amplitudes depend on and, in the case 
of a target polarisation component orthogonal to q, on (fis, the azimuthal angles about 
the virtual-photon direction that are defined in Fig. 2. The dependences on related 
to beam helicity and beam charge have been investigated experimentally [9, 10, 11, 12], 
resulting in first constraints on GPDs. 

At Leading Order (LO) in the fine-structure constant a em , the squared BH 
amplitude |Tbh| is independent of the target polarisation with an unpolarised beam, 
and independent of the lepton charge. In contrast, the squared DVCS amplitude 
|Tdvcs| and the interference term I can depend on the target polarisation even with 
an unpolarised beam, and the sign of the interference term depends also on the lepton 
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Figure 2. Momenta and azimuthal angles for exclusive electroproduction of photons 
in the target rest frame. The quantity <f> denotes the angle between the lepton plane 
containing the three-momenta k and k' of the incoming and outgoing lepton and the 
plane correspondingly defined by q = k — k! and the momentum q' of the real photon. 
The symbol <ps denotes the angle between the lepton plane and S± , the component of 
the target polarisation vector that is orthogonal to q. These definitions are consistent 
with the Trento conventions [13]. 



charge. For an unpolarised lepton beam and a transversely polarised nucleon target, 
these dependences read [14, 15] 
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Here, Sj_ denotes the magnitude of the transverse target polarisation, e\ the beam 
charge in units of the elementary charge, V\ (0)^2(0) contains the ^-dependent lepton 
propagators, and the braces enclose terms that are kinematically suppressed by 1/Q. 
The subscripts 'UU' and 'UT' denote an unpolarised beam with unpolarised and 
transversely polarised targets, respectively. The dependences of the coefficients c n 
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and s n on GPDs are elaborated in Ref. [15], where the kinematic factors K are 
defined. The factor i^ovcs m Eq. 4 suppresses the squared DVCS amplitude by two 
orders of magnitude relative to the interference term in the kinematics of the present 
measurement. Note that the azimuthal angles defined here are different from those used 
in Ref. [15] (0 = 7r — 0[i5] and (f> — (f) s = n + ^[15]), leading to opposite signs for some of 
the coefficients given below. 

The terms of particular interest in this paper appear in bold face in Eqs. 4 and 5. 
The corresponding coefficients can be approximated as 

' f im ins* - en* + isn -n££* 
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with the skewness ^ approximated by C, ~ xb/(2 — %b)- The Compton form factors 
7i, £, 7i and £ are convolutions of hard scattering amplitudes with the corresponding 
twist-two quark GPDs H q , E q , H q and E q , while F\ and F 2 are the nucleon Dirac and 
Pauli form factors [15]. In Eqs. 6-11, the use of bold face differs from that in Eqs. 4 
and 5. Here, the terms not in bold face are suppressed relative to those in bold face in 
the same equation by either xb (or £) or t/M 2 , which are of order 0.1 in the kinematic 
conditions of this measurement. The terms containing xb£ (or ££) are not suppressed 
because the pion-pole contribution to E scales as 1/xb- 

The coefficients c\ vv and c\ vv provide an experimental constraint on the real part 
of the Compton form factors, and can be used to test various models for GPDs as in 
Ref. [12]. Most importantly for the present work, the coefficients c^ut S ' c o,ut an d c i,ut 
provide rare access to the GPD E with no kinematic suppression of its contribution 
relative to those of the other GPDs. Measurements sensitive to these coefficients may 
provide via the Ji relation (Eq. 1) an opportunity to constrain parameterisations of the 
GPD E q in terms of J q [16]. The coefficient c l ovv (cq UT ) has approximately the same 



dependence on GPDs as c\ 



an 



: i,ut) 



The apparent overall suppression of c\ UXJ and 



r 1 

c 0.UT 



by yf—i/Q with respect to c\ VJJ and c\ UT is compensated by an enhancement from 
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y- dependent factors that are not shown, where y = p ■ q/p ■ k. These factors range from 
two to four in the kinematic conditions of the present measurement. The previously 
mentioned strong kinematic suppression of the squared DVCS amplitude relative to the 
interference term is partially compensated by the unshown kinematic factors that apply 
to Eqs. 6, 9 and 10. The net suppression is only about one order of magnitude in the 
HERMES kinematic conditions, and some sensitivity to the GPD E may therefore be 
provided by c^ut S - The coefficient s\ UT provides experimental sensitivity to the GPD 
E, and also to H, which was already probed experimentally through measurements 
of longitudinal target-spin asymmetries [17, 18]. The coefficients cJ?ut S > s °ut S 5 4uui 
c\ UT and S3 UT receive twist-two contributions involving the unknown gluon helicity-ffip 
GPDs [19]. These GPDs do not mix with quark GPDs via Q 2 evolution and thus probe 
the intrinsic gluonic properties of the nucleon [20]. As the contribution of gluon helicity- 
ffip is suppressed by the strong coupling constant as, this contribution competes with 
that from twist-four quark GPDs, which is suppressed by a factor M 2 /Q 2 but not by 
a s [21]. Aside from Cquu S an d C 2UU S ; a ^ other coefficients appearing in Eqs. 4 and 5 are 
related to twist-three quark GPDs. 



3. The experiment 

Hard exclusive production of real photons in the reaction ep^ —>■ e >r yp' is studied. Data 
with a transversely polarised hydrogen target [22] were accumulated using the HERMES 
spectrometer [23] and the longitudinally polarised 27.6 GeV electron and positron beams 
of the HERA accelerator at DESY. This final data set with the transversely polarised 
target was collected over the years 2002-2005. The integrated luminosities for the 
electron and positron samples are approximately 100 pb _1 and 70 pb" 1 , respectively. 

Events are selected if there were detected exactly one photon and one charged 
track identified as the scattered lepton. The hadron contamination in the lepton sample 
is kept below 1% by combining the information from a transition-radiation detector, 
a preshower scintillator detector, and an electromagnetic calorimeter. The kinematic 
requirements imposed are 1 GeV 2 < Q 2 < 10 GeV 2 , 0.03 < x B < 0.35, and v < 22 GeV. 
The real photon is identified through the appearance of a 'neutral signal cluster', which 
is defined as an energy deposition larger than 5 GeV in the calorimeter with a signal 
larger than 1 MeV in the preshower detector, and the absence of a corresponding charged 
track in the back region of the spectrometer. The angular separation # 7 * 7 between the 
virtual and real photons is required to be larger than 5 mrad. This value is determined 
mainly by the lepton momentum resolution. An upper bound of 45 mrad is imposed on 
this angle in order to improve the signal-to-background ratio [24]. 

The recoiling proton was not detected. Instead, an 'exclusive' sample of events 
is selected by requiring the squared missing mass M\ = (q + p — q') 2 of the reaction 
ep — > e'^X to be close to the squared proton mass, where q' is the four-momentum of 
the real photon. The selection criterion was chosen by means of a Monte Carlo (MC) 
simulation of the distribution in M\. The simulation is shown in comparison with 
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Figure 3. Distributions in squared missing-mass from data with positron (filled 
points) and electron (empty circles) beams and from Monte Carlo simulations (solid 
line). The latter include elastic BH (dashed line) and associated BH (filled area) 
processes as well as semi- inclusive background (dotted line). The simulations and data 
are both absolutely normalized. The vertical solid (dashed) lines enclose the selected 
exclusive region for the positron (electron) data. See text for details. 



experimental data in Fig. 3. In the MC simulation [25], the Mo-Tsai formalism [26] 
is used for the elastic BH process that leaves the target nucleon intact, as well as 
the BH process where the nucleon is excited to a resonant state (a category known 
as associated production). For the latter, a parameterisation of the total 7*p cross 
section for the resonance region is used [27]. The individual cross sections for single- 
meson decay channels, e.g., A + — > pit , are treated according to the Maid2000 
model [28]. The remaining contribution is assigned to multi-meson decay channels, 
e.g., A + — > p7r°7T°, whose relative contributions are determined according to isospin 
relations. The simulation also takes into account the semi-inclusive production of neutral 
mesons (mostly tt°) where either only one decay photon is detected or the decay photons 
cannot be resolved. For this, the MC generator Lepto [29] is used in conjunction with 
a set of Jetset [30] fragmentation parameters that had previously been adjusted to 
reproduce multiplicity distributions observed by HERMES [31]. Not shown in Fig. 3 is 
the contribution from exclusive ir° production, which is found to be less than 0.5% in 
the exclusive region using the model in Ref. [32]. The MC yield exceeds the data by 
about 20% in the exclusive region. This may be due to radiative effects not included in 
the simulation, which would move events from the peak to the continuum and improve 
the agreement [33]. On the other hand, if the DVCS process were included in the 
simulation, its contribution would increase the elastic peak. This contribution is highly 
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model-dependent and can vary between 10 and 25% [34]. The exclusive region for the 
positron data is chosen to be —(1.5 GeV) 2 < M\ < (1.7 GeV) 2 , where the value of 
— (1.5 GeV) 2 is displaced from the squared proton mass by three times the resolution in 
M\, and the value of (1.7 GeV) 2 is the point where the contributions from the signal and 
background are equal. As the M\ spectrum of the electron data is found to be shifted by 
approximately 0.18 GeV 2 towards smaller values relative to that of the positron data, 
the exclusive region for the electron data is shifted accordingly. One quarter of the 
effect of this shift on the results presented below is assigned as a systematic uncertainty 
contribution. 

As the recoiling proton remains undetected, t is inferred from the measurement 
of the other final-state particles. For elastic events (leaving the proton intact), the 
kinematic relationship between the energy and direction of the real photon permits t 
to be calculated without using the measured energy of the real photon, which is the 
quantity subject to larger uncertainty. Thus the value of t in the exclusive region is 
calculated as 

1 + h (y - V" 2 + Q 2 cose ri ) ' 

which is exact for elastic events. Using this method, the average resolution (RMS) in t 
is improved from 0.11 to 0.03 GeV 2 . Exclusive events are selected with —t < 0.7 GeV 2 
in order to reduce background. 



4. Azimuthal asymmetries 

Experimental observables that provide sensitivity to the coefficients appearing in Eqs. 4 
and 5 are the Beam-Charge Asymmetry (BCA) 

_ do-+((j>) - da~{<f)) 
Ac((P) = da+(4>) + da-(ct>y (13) 

and the Transverse Target-Spin Asymmetries (TTSAs) 



1 do- + ((f), 4> s ) - da + ((f), 4> s + 7r) + da~((j), (j) S ) - da~((j>, (f) S + 7r) 
S± da+ (0, <j) S ) + da+ (0, 5 + vr) + da~ (<p,<ps) + da~ ((f), 5 + tt) 

1 da + ((p, <f) S ) - da + ((p, 5 + vr) - da~((p, <f> s ) + da~((f), <p s + tt) 



(14) 



(15) 



S± da + ((f),(j)s) + da + ((j),(f) S + n) + da (0, (j> s ) + da (0, (f) S + n) 
Here the subscripts on the A's represent dependence on beam Charge (C) or Transverse 
(T) target polarisation, with an Unpolarised (U) beam, and the superscripts ± stand for 
the lepton beam charge. These asymmetries are related to the coefficients in Eqs. 3-5 
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by: 
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HERMES reported the first measurement of the BCA [12], providing access to c* UTJ , 
the coefficient of the cos0 modulation of the interference term for an unpolarised target. 
The present paper reports more precise BCA results using a considerably larger new 
data set from the first DVCS measurement done with transverse target polarisation. 
Most importantly, the extracted TTSAs provide access to Cqut S > c out an d c\ UT , which 
are sensitive to the total angular momentum of quarks in the nucleon, as noted above. 



4-1. Extraction of azimuthal asymmetry amplitudes 
The distribution in the expectation value of the yield is given by: 

(N)(S±, e h 0, S ) = £ (S±, ei) e(e h 0, 5 ) a vu {4>) x 

+ 5 ± ^ cs (0, S ) + e^ c (0) + ei S ± A\j 



(22) 



Here £ is the integrated luminosity, e the detection efficiency, and auu the cross section 
for an unpolarised target averaged over both beam charges. The BCA «4c(0) and the 
TTSAs ^4u^ cs (0,0s) and *4{j T (0, S ) in Eq. 22 are expanded in terms of the same 
harmonics in and 0s as those appearing in Eqs. 4 and 5 (as well as the harmonics 
cos(0 — 0s) sin(30) and sin(0 — 0s) cos(30) in Eq. 24, included as a systematic check): 
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^ cs (0, 05; ^ cs ) = E ^ T ( y c t os( ^ sin(0 - S ) cos(n0) 

n=0 

+ E ^utDVCS )sin(n0) - S ) sin(n0), (24) 



n=l 



3 



4jt(0, 05; ffc) = E ^-^ )cos( ^ sin(0 - S ) cos(n0) 

n=0 

3 

+ A™ { *-* s)sia{n<i,) cos(0 - S ) sin(n0). (25) 

n=l 

Here 77c, T7ut CS an< ^ ^ut represent the sets of Fourier coefficients or azimuthal 
asymmetry amplitudes, hereafter called 'asymmetry amplitudes', appearing in the right- 
hand sides of Eqs. 23-25 describing respectively the dependences of the squared DVCS 
amplitude and interference term on beam Charge (C), Transverse (T) target polarisation 
or both, with an Unpolarised (U) beam. These 18 asymmetry amplitudes embody the 
essential sensitivities to GPD models of the coefficients of the corresponding functions 
of appearing in Eqs. 4 and 5, to the degree that one can neglect the effects of the 
coefficients c^uu an d c 2 uu anc ^ the squared unpolarised DVCS amplitude in Eqs. 18-20 
and the effect of the 0-dependence of the BH propagators. In any case, the extracted 
asymmetry amplitudes are well defined and can be computed in various GPD models 
for direct comparison with the data. For each kinematic bin in —t, xb or Q 2 , they are 
simultaneously extracted from the observed exclusive event sample using the method of 
Maximum Likelihood. The distribution of events is parameterised by the function Af par , 
which is defined as 

■A/" P ar(<Sj_, e h 0, (f) S ; ?7ut CS > Vc, "Hxjt) = £ ( s ±, e i) e ( e h 0, <j>s) vuu(4>) x 

1 + S ± A^ CS (<P, 5 ; r?OT CS ) + eiAcifa rjc) + e^ ± ^ T (0, 5 ; »£ T )] • (26) 

While the net beam polarisations of both positron and electron data samples used in 
the current measurement are not completely negligible (0.03 ± 0.02 and —0.03 ± 0.02, 
respectively), algebraic investigations (confirmed by MC studies) show that this does 
not affect the asymmetry amplitudes presented here. Not included in Eqs. 23-25 are 
negligible terms involving the small component of the target polarisation that is parallel 
to q [35]. 

Within the scheme known as Extended Maximum Likelihood [36], the likelihood 
function L to be minimized is taken as 

- In L(/7^ cs , rye, T7ut) = Afpar^uT^ Vc, vIjt) 

- E In [l + S^ cs ^\ 5 ; 77^ cs ) + e^ c (0 4 ; Vc) 

i=l 

+ e^K T (0\0 5 ; 77^)1, (27) 
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where iV is the observed number of events, and the parameterised total number A/" par 
of events is 

Ar p ar(r7u¥ CS ,r7c,r7u T ) = 

I dS ± d<j> d<p s KAS±, e h 0, 4> s] 77^ cs , rye, r}]j T ). (28) 
J ei =±i 

The cross section auu and the detection efficiency e do not depend on rj^ cs , r] c or 
T7[j T and thus cannot affect the location of the likelihood maximum. Hence they have 
been omitted in the logarithms in Eq. 27. It is also not necessary to consider them 
explicitly in evaluating A / par (r7u^ cs , r] c , r]\j T ) in Eq. 28, because the needed information 
about them is encoded in the total yields obtained by combining events for both 
target polarisations and beam charges. Luminosity imbalances between beam charges 
or target polarisations are taken into account by assigning weights Wi to the events, 
which are adjusted to provide effectively vanishing net target polarisation and net beam 
charge for this combined data set. The weights are normalized to also retain the same 
integrated luminosity Ctot as the observed 'unweighted' data sample. The resulting 
event distribution corresponds to the product C t at e (<P, ( f ) s) o'uui.'P)- m this manner, an 
event distribution corresponding to Eq. 26 is constructed to estimate the parameterised 
total number of events in Eq. 28: 

i=i 



4-2. Background corrections and systematic uncertainties 

The results from the minimization of Eq. 27 in each kinematic bin are corrected for 
background from semi-inclusive and exclusive production of neutral mesons, mainly 
pions, in order to estimate the true asymmetry amplitude: 

A t = ^— -, 29 

1 — s — e 

where A r stands for the extracted raw asymmetry amplitude, and s and A s (e and 
A e ) the fractional contribution and corresponding asymmetry amplitude of the semi- 
inclusive (exclusive) background. The combination of these background contributions 
s + e ranges from 2 ± 1% to 11 ± 5% in the kinematic space [34]. As these background 
contributions are only very weakly beam-charge dependent, their asymmetries with 
respect to the beam charge or to the product of the beam charge and the transverse 
target polarisation are neglected. The asymmetry of the semi-inclusive ir° background 
with respect to only the transverse target polarisation is extracted from experimental 
data by requiring two 'neutral signal clusters' in the calorimeter with their invariant mass 
between 0.10 and 0.17 GeV. The restriction on the energy deposition in the calorimeter 
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of the less energetic neutral signal cluster is relaxed to 1 GeV to improve the statistical 
precision. The fractional energy of the reconstructed neutral pions is required to be 
large, z = E n /u > 0.8, as only these contribute to the exclusive region according 
to MC simulations [34]. These simulations showed that the extracted 7r° asymmetry 
does not depend on whether only one or both photons are in the acceptance. It is 
convenient to use the direction of the reconstructed pion in place of that of the photon 
to calculate the azimuthal angles <fi and For the exclusive 7r° background, the 
asymmetry amplitudes with respect to only target polarisation are not extracted due 
to the limited statistical precision but rather assumed to be ± 1. After applying 
Eq. 29, the resulting asymmetry amplitude A t is expected to originate from only elastic 
and associated production. On average 12% of the BH cross section arises from the 
latter [34], according to the simulation described above. The kinematic dependences of 
this contribution are shown in Fig. 4. No correction is made or uncertainty assigned for 
associated production, as it is considered to be part of the signal. 




1 io 



0.2 0.4 0.6 0.1 0.2 0.3 2 4 6 8 10 



■t (GeV z ) 



Q 2 (GeV 2 ) 



Figure 4. Kinematic dependences of the simulated fractional contributions from 
associated production. See text for details. 



The dominant contributions to the total systematic uncertainty are those from 
the detector acceptance and finite bin width, and the determination of the target 
polarisation [34]. The combined contribution to the systematic uncertainty from the 
detector acceptance, finite bin width, and the alignment of the detector elements with 
respect to the beam, including possible effects from the beam and track curvature in 
the transverse magnetic field of the target, is determined from MC simulations based on 
five GPD models described in Ref. [37]. In each kinematic bin, it is defined as the RMS 
difference between the asymmetry amplitude extracted from the MC data in that bin by 
minimizing Eq. 27 and the corresponding model predictions calculated analytically at 
the mean kinematic values of that bin given in Table 1. The other sources are associated 
with the background correction, calorimeter calibration and the relative shift of the M\ 
spectra between the positron and electron data. These contributions, given in Table 2, 
are added in quadrature to form the total systematic uncertainty per kinematic bin, 
appearing in Table 1. Not included is any contribution due to additional QED vertices, 
as the most significant of these has been estimated to be negligible [38]. 
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Figure 5. Asymmetry amplitudes describing the dependence of the interference term 
on the beam charge (^4c)j for the exclusive sample. The squares represent the results 
from the present work, while data represented by triangles (shifted right for visibility) 
were reported in Rcf. [12]. The filled symbols indicate those results of greatest interest 
(see text). The error bars (bands) represent the statistical (systematic) uncertainties. 
The curves are predictions of variants of a double-distribution GPD model [32, 39], 
with profile parameters given in Table 3. See text for details. 



5. Results 

Figures 5-7 show as a function of —t, Xb or Q 2 , in four bins, the results from the 
combined fit. The 'overall' results in the left-most columns correspond to the entire 
experimental acceptance. Fig. 5 shows the amplitudes related to only beam charge, 
while Figs. 6 and 7 show the amplitude Aut,dvcs, which relates to transverse target 
polarisation only, and Aut\i, which relates to both. For simplicity of presentation, the 
amplitudes ^4ut,dvcs an d ^4ut,i fc> r the same azimuthal dependence are shown together 
in each panel, even though they typically do not relate to the same GPDs. The filled 
symbols represent the asymmetry amplitudes of interest here (see Table 1), related to 
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Figure 6. Asymmetry amplitudes describing the dependence of the squared DVCS 
amplitude (circles, Aut.dvcs) and the interference term (squares, Atjt.i) on the 
transverse target polarisation, for the exclusive sample. The filled symbols indicate 
those results of greatest interest (see text). The circles (squares) are shifted right 
(left) for visibility. The error bars represent the statistical uncertainties, while the top 
(bottom) bands denote the systematic uncertainties for Atjt.i (^.ut,dvcs)> excluding 
the 8.1 % scale uncertainty from the target polarisation measurement. The curves are 
predictions of the GPD model variant (Reg, no D) shown in Fig. 5 as a continuous 
curve, with three different values for the u-quark total angular momentum J u and fixed 
(i-quark total angular momentum Jd = [16]. See text for details. 



the coefficients given in Eqs. 6-11, of the corresponding harmonics of appearing in 
Eqs. 4 and 5. 

Of particular interest is the asymmetry amplitude A c ^ s ^ in the upper row of Fig. 5. 
Equation 7 shows that this amplitude is sensitive to the GPD H in the HERMES 
kinematic conditions. Also shown in this figure is the previously published result, which 
has been shown to constrain GPD models [12]. The greatly improved precision of 
the present measurement confirms that this amplitude increases with increasing —t. As 
mentioned above regarding the corresponding coefficients c\ UTJ and c\ ULr , the amplitude 
^coso</> j g eX p ec t ec } to relate to the same combination of GPDs as does A™ s *. The results 
shown in Fig. 5 suggest that the magnitude of this amplitude also increases with —t, 
while its opposite sign is expected from Eq. 8. 

Of special interest in this work are the amplitudes ^uti ^ s ^ cos ^\ n = 0,1, 
presented in the top two rows of Fig. 6. Equations 10 and 9 show that these amplitudes 
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Figure 7. Asymmetry amplitudes that are expected to be suppressed, presented as 
in Fig. 6, except that the curves are calculated only for J u = 0.4. 



are sensitive to the GPD E and hence to the total angular momenta of quarks. These 
amplitudes are found to have substantial magnitudes with opposite signs but little 
kinematic dependence, possibly increasing in magnitude with —t. Their opposite signs 
are expected from Eq. 10. Also of interest is the amplitude ^utdvcS' shown in the top 
row of Fig. 6, which Eq. 6 suggests is also sensitive to the GPD E. The overall result is 
non-zero by 2.8 times the total uncertainty. These data tend to increase in magnitude 
at larger values of Q 2 . (In fixed-target experiments, xb and Q 2 are strongly correlated.) 
They provide the first experimental evidence for an azimuthal harmonic in the squared 
DVCS amplitude, in this case related to transverse target polarisation. 

The amplitude ^^(^"^sm^ Su0wn j n bottom row of Fig. 6 is sensitive mainly 

to the GPDs H and E, while the contribution from the GPD E is suppressed by an 
additional factor of xb (see Eq. 11). The measured asymmetry amplitudes are consistent 
with zero. 

The amplitudes represented by the unfilled symbols are expected to be suppressed, 
and are indeed found to be typically small. However, values that depart from zero by 
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more than twice the total uncertainty are found for the entire experimental acceptance 
for two of the four amplitudes in Fig. 7 that receive a contribution from gluon helicity- 
flip, which are A'^^vcs an< ^ ^ut i ^ s ^ sm ^ . The asymmetry amplitudes related 

to the squared DVCS amplitude in the bottom two rows of Fig. 7 correspond to 
coefficients that do not appear in Eq. 4 as a consequence of the one-photon exchange 
approximation. They are found to be consistent with zero. 

6. Comparison with theory 

The data are compared with various theoretical calculations to LO in a em and a s , 
which do not account for the contributions of associated production. The calculations 
shown in Figs. 5-7 employ variants of a GPD model developed in Refs. [32, 39]. 
These are based on the widely used framework of double distributions involving a 
product of PDFs representing the forward limit and a profile function representing the 
skewness dependence [40]. The forward limits of the GPDs H q are constructed using 
the MRST98 [41] parameterisation of PDFs evaluated at the measured Q 2 value for 
each data point. More modern parameterisations are expected to result in a negligible 
difference [16]. In the model description for H and E, the forward limit of the GPD H 
is fixed by the quark helicity distributions Aq(xs,Q 2 ), while the GPD E is evaluated 
from the pion pole, which provides only the real part. The 'profile parameters' b va i 
and b sea control the skewness dependence of GPDs for the valence and sea quarks, 
respectively [40, 42]. The t dependence of the GPDs is calculated in either the simplest 
ansatz where the t dependence factorises from the t-independent part H q (x,£), or in 
the Regge-inspired ansatz. The GPDs H and E are optionally modified by the so-called 
(t-independent) D-term [43] contribution to the double-distribution part of the GPD, 
with a value calculated in the chiral-quark soliton model [44]. The twist-three GPDs 
are treated in the Wandzura-Wilczek approximation, and the gluon helicity-flip GPDs 
are not included. The quark total angular momenta J q of quarks and antiquarks of 
flavour q (q = u,d) enter as model parameters for the GPD E. The strange sea is 
neglected. The computational program of Ref. [45] is used. The calculation is done 
at the average kinematics of every bin (see Table 1). For the comparison of the BCA 
amplitudes to the double-distribution model shown in Fig. 5, the model variations of 
interest are those that change the GPD H, since the impact of the GPDs H and E is 
suppressed at HERMES kinematic conditions (see Eq. 7). Four different variants are 
selected by choosing either a factorised (Fac) or a Regge-inspired (Reg) t dependence, 
each with or without the contribution of the D-term. While these four variants lead to 
very different model predictions for A C Q Sr ^ as illustrated in the figure, the variation of 
the profile parameters b va i and b sea lead to smaller changes [24]. However, by comparing 
the data for the cos <fi amplitude with predictions of all four variants of this model in 
combination with four specific sets of values for the profile parameters, it is found that 
the calculation using the Regge-inspired t dependence without the D-term and b va i = oo, 
b sea = 1 results in a confidence level much higher than all the alternatives. Here, b = 1 
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corresponds to a substantial skewness dependence, which is eliminated for b = oo. All 
the variants shown in Fig. 5 are calculated using profile parameter values that yield 
the best agreement with data (see Table 3). For any choice of the profile parameters, 
the variant Regge with a D-term is excluded, while the factorised ansatz is disfavoured 
either with or without the D-term. The factorised t dependence is also disfavoured on 
theoretical grounds [46, 47]. 

The theoretical calculations of the TTSA amplitudes shown in Figs. 6 and 7 are 
made using the Regge inspired t dependence without the D-term and b va i = oo and 
bsea = 1, a combination that is favoured by the BCA data as described above. However, 
while the calculated TTSA amplitudes are less sensitive to that choice, some of them 
are quite sensitive to the choice of the quark total angular momenta J q [16]. This 
sensitivity is illustrated by the curves in Fig. 6 evaluated with three different values 
of J u (0.2, 0.4, 0.6), while fixing Jd = [16]. Although this model fails to describe 
the data for Ay-p ^y^g j 

the model curves confirm the expectation from Eqs. 6, 9 and 
10 that the TTSA amplitudes ^t,dvcs; Af^ s)cosrp and Af^ s) have significant 
sensitivity to J u . However, for this double-distribution model, the amplitudes related to 
the interference term are in reasonable agreement with the data, of which A s ^\ t t's)c°s<t> 
has the greatest sensitivity. The curves in Fig. 5 and 7 are evaluated with fixed J u = 0.4 
and Jd = 0, since the sensitivity here to J u and Jd is negligible. 

The BCA amplitude A c ° s * and the TTSA amplitudes Af^ s) and A s j^ s)coB0 
are also compared to calculations based on the 'dual-parameterisation' model of 
GPDs [48, 49] in Figs. 8 and 9. Calculations for all other amplitudes are not shown 
since the model contains neither GPDs H and E nor higher-twist contributions. For 
the BCA amplitude, the curves in Fig. 8 are evaluated with fixed J u = Jd = 0, 
which best describes the TTSA amplitudes as discussed below. The t dependence 
of the GPDs is assumed to be either factorised and exponential, or non-factorised 
and Regge-inspired [39]. Both choices describe the BCA data equally well, but with 
a smaller confidence level (x 2 /d.o.f.=2.2 for both Fac and Reg) than the favoured 
double-distribution model described above. On the other hand, the existing beam-spin 
asymmetry data [51] are better described by this dual-parameterisation model. The 
Regge-inspired variant of this model with J u = 0, 0.2, and 0.4, and Jd = is used for 
the calculations shown in Fig. 9. It is apparent that also in the dual-parameterisation 
model, the amplitudes ^^-"fe anc j A s ^\~^ are sensitive to J u . 

7. Quark total angular momentum 

In either GPD model discussed above, the GPD E is parameterised in terms of J u 
and Jd [39, 16]. In both cases these two parameters are fit to the measured overall 
TTSA amplitudes A™ ( f^ )cos0 and appearing in the left column of Figs. 6 

and 9, respectively. The other parameters for the respective GPD models are the same 
as used for the curves in Figs. 6 or 9. The area in the (J u , Jd)-plane in which the 
reduced X 2 ~X m in value is not larger than unity corresponds to a one- standard- deviation 
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Figure 8. Similar to the top row of Fig. 5, except that the curves are calculations 
based on the dual-paramctcrisation GPD model [49]. See text for details. 
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Figure 9. Similar to top two rows of Fig. 6, except that the curves are the 
calculations [50] using the Regge-inspired form of the t dependence in the dual- 
paramcterisation GPD model [49] (shown in Fig. 8 as a dash-dotted curve) . with three 
different values for the w-quark total angular momentum J u and fixed d-quark total 
angular momentum Jd = 0. See text for details. 



constraint on J u vs Jd- This area is found to be one of the sloped bands in Fig. 10, 
which in units of h can be represented in the case of the double-distribution model as 

J u + J d /2.8 = 0.49 ± 0.17(exp tot ), (30) 

and in the case of the dual-parameterisation model as 

J u + J d /2.8 = -0.02 ± 0.27(ex Ptot ). (31) 

The uncertainty is propagated from the total experimental uncertainty in the measured 
TTSA amplitudes. This uncertainty dominates the effects of variations within either 
of the GPD models of the values of profile parameters b and the inclusion of the D 
term. The variation of the value of J u + Jd/2.8 obtained by fitting the amplitudes 
in the bins of the three kinematic variables is found to be not larger than ±0.15. The 
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Figure 10. Model-dependent constraints on u-quark total angular momentum J u 
vs d-quark total angular momentum Jd, obtained by comparing DVCS experimental 
results and theoretical calculations. The constraints based on the HERMES data 
for the TTSA amplitudes 4jM*-*0«»* &nd A *Ht-<t>s) use the double-distribution 
(HERMES DD) [32, 39] or dual-parametcrisation (HERMES Dual) [49] GPD models. 
The additional band (JLab DD) is derived from the comparison of the double- 
distribution GPD model with neutron cross section data [55]. Also shown as small 
(overlapping) rectangles arc results from lattice gauge theory by the QCDSF [52] and 
LHPC [47] collaborations, as well as a result for only the valence quark contribution 
(DFJK) based on zero-skewness GPDs extracted from nucleon form factor data [53, 54]. 
The sizes of the small rectangles represent the statistical uncertainties of the lattice 
gauge results, and the parameter range for which a good DFJK fit to the nucleon form 
factor data was achieved. Theoretical uncertainties are unavailable. 



results from fitting the two azimuthal amplitudes separately were found to be consistent. 
The large difference between the constraints obtained using the double-distribution 
and dual-parameterisation models is an indication of a large model dependence of the 
{Jui Jd) constraint obtained, which may be related to the failure of both models to fully 
describe all other available DVCS data. Both constraints are consistent with results, 
also shown in Fig. 10, from unquenched lattice gauge simulations by the QCDSF [52] 
and the LHPC [47] collaborations. The statistical uncertainties of the lattice gauge 
results are comparable to the size of the plotted symbols. The QCDSF calculation of 
the first moments of the GPDs, the so-called generalized form factors, is based on a 
simulation using dynamical Wilson fermions with pion masses down to 350 MeV. The 
dynamical LHPC calculation is based on a hybrid approach of rooted staggered sea and 
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domain wall valence quarks. In both calculations the generalized form factors have been 
simultaneously extrapolated in t and to t = and the physical point, respectively, 
using the same results from chiral perturbation theory. Both calculations include only 
contributions from connected diagrams. The uncertainties are primarily statistical but 
include some systematic uncertainties from the chiral, continuum and infinite volume 
extrapolations. Also shown in Fig. 10 is a result for only the valence contribution 
to the quark total angular momenta [53, 54]. It is based on the extraction of zero- 
skewness GPDs from nucleon form factor data, assuming handbag diagram dominance 
and exploiting well known sum rules. The size of the plotted symbol corresponds to the 
parameter range for which a good fit to the nucleon form factor data was achieved. 



Table 1. Results of particular interest for the asymmetry amplitudes of the 
asymmetries with respect to the beam charge and transverse target polarisation for 
the exclusive sample. 
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-0.105 ±0.047 ±0.013 
-0.201 ±0.058 ±0.027 


> 

CD 

o 
Q? 


1.0-1.5 
1.5-2.3 
2.3-3.5 
3.5-10.0 


0.08 
0.10 
0.13 
0.19 


0.06 
0.08 
0.11 
0.17 


1.2 
1.9 
2.8 
4.9 


-0.014 ±0.019 ±0.016 
-0.004 ±0.018 ±0.016 
-0.023 ±0.021 ±0.015 
-0.003 ±0.023 ±0.016 


0.025 ±0.026 ±0.011 
0.070 ± 0.026 ± 0.015 
0.058 ± 0.030 ± 0.008 
0.005 ± 0.032 ± 0.014 


0.044 ±0.056 ±0.012 
-0.080 ±0.046 ±0.010 
-0.113 ±0.049 ±0.012 
-0.143 ±0.048 ±0.015 



kinematic bin 


(-*) 
(GeV 2 ) 


(x B ) 


(Q 2 ) 

(GeV 2 ) 


.sin (0-0 s ) 
A UT, I 

±<5stat ± <5 sys t 


.sin (0 — S ) cos 
A UT, I 

±<5stat ± <5syst 


. cos (0 — S ) sin 
rt UT, I 

±5stat ± <5 s yst 




overall 


0.12 


0.09 


2.5 


0.035 ± 0.024 ± 0.024 


-0.164 ±0.039 ±0.023 


0.005 ±0.040 ±0.015 


CI 

> 

CD 

O 
-+~ 

i 


0.00-0.06 
0.06-0.14 
0.14-0.30 
0.30-0.70 


0.03 
0.10 
0.20 
0.42 


0.08 
0.10 
0.11 
0.12 


1.9 
2.5 
2.9 
3.5 


-0.030 ±0.031 ±0.008 
0.022 ±0.044 ±0.021 
0.133 ±0.050 ±0.025 
0.085 ± 0.082 ± 0.028 


-0.152 ±0.068 ±0.026 
-0.073 ± 0.068 ± 0.008 
-0.244 ± 0.078 ± 0.028 
-0.294 ±0.126 ±0.026 


-0.100 ±0.069 ±0.044 
0.054 ± 0.076 ± 0.030 
0.144 ±0.083 ±0.020 
0.024 ±0.113 ±0.029 


H 


0.03-0.07 
0.07-0.10 
0.10-0.15 
0.15-0.35 


0.10 
0.10 
0.13 
0.20 


0.05 
0.08 
0.12 
0.20 


1.5 
2.2 
3.1 
5.0 


0.083 ±0.051 ±0.021 
0.037 ± 0.048 ± 0.021 
-0.033 ±0.048 ±0.021 
0.048 ± 0.055 ± 0.024 


-0.166 ±0.084 ±0.047 
-0.148 ±0.078 ±0.034 
-0.100 ±0.072 ±0.020 
-0.182 ±0.084 ±0.026 


-0.034 ±0.081 ±0.025 
-0.078 ±0.080 ±0.015 
0.078 ± 0.073 ± 0.025 
0.066 ±0.088 ±0.056 


> 

CD 

o 


1.0-1.5 
1.5-2.3 
2.3-3.5 
3.5-10.0 


0.08 
0.10 
0.13 
0.19 


0.06 
0.08 
0.11 
0.17 


1.2 
1.9 
2.8 
4.9 


0.117 ±0.056 ±0.024 
-0.043 ± 0.046 ± 0.026 

0.066 ± 0.049 ± 0.028 
-0.002 ± 0.049 ± 0.020 


-0.174 ±0.092 ±0.047 
-0.170 ±0.073 ±0.031 
-0.249 ±0.078 ±0.025 
-0.059 ±0.072 ±0.011 


-0.034 ±0.093 ±0.018 
-0.036 ± 0.079 ± 0.020 
0.028 ± 0.076 ± 0.026 
0.056 ± 0.079 ± 0.035 
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8. Conclusions 

Transverse target-spin azimuthal asymmetries in electro-production of real photons 
are measured for the first time, and for both beam charges. A combined fit of 
this data set separates for the first time the azimuthal harmonics of the squared 
DVCS amplitude and the interference term. The extracted charge asymmetry of 
the interference term is much more precise than previously published results, and 
constrains models for Generalized Parton Distributions. By comparing GPD-model 
calculations with extracted azimuthal asymmetry amplitudes associated with both 
beam charge and transverse-target polarisation, a model-dependent constraint on 
the total angular momenta carried by u and <i-quarks in the nucleon is obtained 
as J u + Jd/2.8 = 0.49 ± 0.17(exp tot ) using a double-distribution GPD model, and 
Ju + Jd/2.8 = —0.02 ± 0.27(exp tot ) using the dual-parameterisation model. Thus, such 
data have the potential to provide quantitative information about the spin content of 
the nucleon when GPD models become available that fully describe all existing DVCS 
data. 

Table 2. Systematic uncertainties of the results of particular interest for the azimuthal 
amplitudes of the asymmetries with respect to the beam charge and transverse 
target polarisation for the exclusive sample. Those results involving transverse 
target polarisation are also subject to an additional 8.1% scale uncertainty from the 
determination of the target polarisation. 
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M% shift 


0.004 


0.001 


0.000 


0.000 


0.001 


0.001 


Background correction 


0.000 


0.001 


0.005 


0.001 


0.004 


0.000 


Calorimeter calibration 


0.001 


0.003 


0.002 


0.002 


0.005 


0.003 


Acceptance, bin width, alignment 


0.017 


0.015 


0.002 


0.024 


0.019 


0.014 
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Table 3. The quality of agreement between the measured A c ^ ' asymmetry 
amplitude and four variants of the double-distribution model. These \ 2 values are 
based on sums over the 12 bins in —t, xb and Q 2 , without considering the correlations 
among them because of their sharing of some events. The values given are the sums 
divided by 12. Values shown in bold face denote for each variant the profile parameter 
set that yields the best agreement with the data, chosen to produce the curves in Fig. 5. 



X 2 /d.o.f. (J u = 0.4) 


b v = b s = 1 
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b v = oo, b s = 1 


by = b s = oo 
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1.1 
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